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Abstract 

A new class of a spatially homogeneous and anisotropic Bianchi type-I cosmological models of the universe 
for perfect fluid distribution within the framework of scalar-tensor theory of gravitation proposed by Saez and 
Ballester (Phys. Lett. 113:467, 1986) is investigated. To prevail the deterministic solutions we choose the 
different scale factors which yield time-dependent deceleration parameters (DP) representing models which 
generate a transition of the universe from the early decelerated phase to the recent accelerating phase. Three 
different physically viable models of the universe are obtained in which their anisotropic solutions may enter 
to some isotropic inflationary era. The modified Einstein's field equations are solved exactly and the models 
are found to be in good concordance with recent observations. Some physical and geometric properties of 
the models are also discussed. 

Key words: Bianchi type-I universe. Exact solution. Alternative gravitation theory. Variable deceleration pa- 
rameter 
PACS: 98.80.-k 

1 Introduction 

General relativity (GR) passes all present tests with flying colours, however, there are several reasons why it 
remains very important to consider alternative theories of gravity. The first one is that theoretical attempts 
at quantizing gravity or unifying it with other interactions generically predict deviation from Einstein's theory, 
because gravitation is no longer mediated by a pure spin-2 field but also by partners to the graviton. The second 
reason is that it is any way extremely instructive to contrast GR's predictions with those of alternative models, 
even if there were no serious theoretical motivation for them. The third reason is the existence of several puzzling 
experimental issues, which do not contradict GR in a direct way, but may nevertheless suggest that gravity does 
not behave at large distances exactly as Newton and Einstein predicted. Cosmological observations notably tell 
us that about 96% of the total energy density of the universe is composed of unknown, non-baryonic, fiuids (72% 
of "dark energy" and 24% of "dark matter" ) , and the acceleration of the two Pioneer spacecrafts towards the 
Sun happens to be larger than what is expected from the law. Therefore, alternative theories are any way 
important to study. 



Since the observed universe is almost homogeneous and isotropic, space-time is usually described by a 
Friedman-Lemaitre-Robertson- Walker (FLRW) cosmology. But it is also believed that in the early universe 
the FLRW model does not give a correct matter description. The anomalies found in the cosmic microwave 
background (CMB) and the large structure observations stimulated a growing interest in anisotropic cosmo- 
logical model of the universe. Observations by the Differential Macrowave Radiometers (DMR) on NASA's 
Cosmic Background Explorer Spacecraft registered anisotropy in various angle scales. It is conjectured, that 
these anisotropies hide in their hearts the entire history of the cosmic evolution down to recombination, and they 
are considered to be indicative of the universe geometry and the matter composing the universe. It is expected, 
that much more will be known about anisotropy of cosmic microwave's background after the investigations of 
the microwave's anisotropy probe. There is a general agreement among cosmologists that cosmic microwave's 
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background anisotropy in the small angle scale holds the key to the formation of the discrete structure. The 
theoretical argument [T] and the modern experimental data support the existence of an anisotropic phase, which 
turns into an isotropic one. 

In recent years, there has been considerable interest in scalar tensor theories of gravitation which are con- 
sidered to be essential to describe the gravitational interaction near the planks scale, string theory, extended 
inflation and many higher order theories imply scalar field. Scalar- Tensor theories of gravitation provide the 
natural generalizations of general relativity and also provide a convenient set of representations for the obser- 
vational limits on possible deviations from general relativity. There are two categories of gravitational theories 
involving a classical scalar field (f>. In first category the scalar field (p has the dimension of the inverse of the 
gravitational constant G among which the Brans-Decke theory [5] is of considerable importance and the role of 
the scalar field is confined to its effect on gravitational field equations. Brans and Decke formulated a scalar- 
tensor theory of gravitation which introduces an additional scalar field besides the metric tensor gij and a 
dimensionless coupling constant uj. This theory goes to general relativity for large values of the coupling constant 
uj > 500. In the second category of theories involve a dimensionless scalar field. Sdez and Ballester [5] developed 
a scalar-tensor theory in which the metric is coupled with a dimensionless scalar field in a simple manner. This 
coupling gives a satisfactory description of the weak fields. In spite of the dimensionless character of the scalar 
field, an anti-gravity regime appears. This theory suggests a possible way to solve the missing-matter problem 
in non-flat FRW cosmologies. The Scalar- Tensor theories of gravitation play an important role to remove the 
graceful exit problem in the inflation era In earlier literature, cosmological models within the framework 
of Saez-Ballester scalar-tensor theory of gravitation, have been studied by Singh and Agrawal [3 E] , Ram and 
Tiwari [7] , Singh and Ram [8 . Mohanty and Sahu |9l [10] have studied Bianchi type-VIo and Bianchi typc-I 
models in Saez-Ballester theory. In recent years, Tripathi et al. [TT], Reddy et al. [T21II3], Reddy and Naidu ^14j . 
Rao et al. [ll[T6l[n], Adhav et al. [H], Katore et al. [ig, Sahu [20], Singh [21] and Pradhan and Singh [22] 
have obtained the solutions in Sdez-Ballester scalar-tensor theory of gravitation in different context. Recently, 
Naidu et al. [23] and Reddy et al. [24] have studied LRS Bianchi type-II and five dimensional dark energy 
models in Sdez and Ballester scalar tensor theory of gravitation respectively. 

Recently, Kumar and Singh ,,25) obtained exact Bianchi type-I cosmological models in Sdez and Ballester 
Scalar- Tensor theory of gravitation by assuming the constant deceleration parameter. In literature it is common 
to use a constant deceleration parameter, as it duly gives a power law for metric function or corresponding quan- 
tity. Several cosmological observations indicate that the observable universe is undergoing a phase of accelerated 
expansion (Riess et al. [2^; Perlmutter et al. [27]; Bennet et al. [28]; Tegmark et al. [29]; Allen et al. [SOJ. 
Recent observations of SNe la of high confidence level (Tonry et al. 01]; Riess et al. [31]; Clocchiatti et al. [55] ) 
have further confirmed this. Also, the transition redshift from deceleration expansion to accelerated expansion 
is about 0.5. Now for a Universe which was decelerating in past and accelerating at the present time, the DP 
must show signature flipping (see the Refs. Padmanabhan and Roychowdhury [34], Amendola [35], Riess et al. 
[36]). So, in general, the DP is not a constant but time variable. 

Socorro et al. [37] and Jamil et al. [38" have studied on anisotropic Bianchi type-I cosmology in Sdez-Ballester 
scalar-tensor theory of gravity. Recently, Pradhan et al. [39] have studied some exact Bianchi type-I cosmo- 
logical models in scalar-tensor theory of gravitation with time dependent deceleration parameter. Motivated by 
the discussions, in this paper, we have obtained a new class of exact solutions of field equations given by Sdez 
and Ballester |3j in Bianchi type-I space-time by considering a time dependent deceleration parameter. The out 
line of the paper is as follows: In Section 2, the metric with a Lagrangian and the field equations are described. 
Section 3 deals with the solutions of the field equations. Subsections 3.1, 3.2 and 3.3 describe the three different 
models by considering the three different values of scale factors which yield the time dependent deceleration 
parameters and physical and geometric behaviour of the model are also discussed. Finally, conclusions are sum- 
marized in the last Section 4. 



2 The Metric and Field Equations 

We consider spatially homogeneous and anisotropic Bianchi type-/ space-time given by 

ds^ ^ -dt^ + A^dx^ + B'^dy'^ + C'^dz^, (1) 
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where the metric potentials A, B and C are functions of cosmic time t alone. This ensures that the model is 
spatially homogeneous. 



We, consider the simple case of a homogeneous but anisotropic Bianchi type-I model with matter term with 
a scalar field </). Our model is based on a non-standard scalar-tensor theory, defined in Sdez and Ballester [3] 
with a dimensionless scalar field </> and tensor field gij . This alternative theory of gravitation is combined scalar 
and tensor fields in which the metric is coupled with a dimensionless scalar field. We assume the Lagrangian 

i = i? - (2) 

R being the scalar curvature, <j) a dimensionless scalar field, uj and k arbitrary dimensionless constants and 
the contraction (p^ag"'^ ■ Here the partial derivatives are denoted by a comma and the covariant derivatives by a 
semicolon in the usual way. 

From the above Lagrangian we can establish the action 

/ = / {L + xLm){-g)^dx^dx^dx^dx'^, (3) 

where L„i is the matter Lagrangian, g is the determinant of the matrix the coordinates, 

an arbitrary region of integration and x = Stt (we use geometrized units). When fc = 0, our model is just the 
Einstein gravity with a massless minimally coupled scalar field coupled to gravity. Different aspects of this model 
have been investigated in the literatures (already described in Sect. 1). By considering arbitrary independent 
variations of the metric and the scalar field vanishing at the boundary of the variation principle 

51 = 0, (4) 

leads to a generalized Einstein equation 

2p''^:; + kp''-^<f>,iP'' ^0, (5) 

where Gij = Rij — \Rgij is the Einstein tensor; Ty is the stress-energy tensor of the matter Lagrangian Lm. 



Since the action / is a scalar, it can be easily proved that the equation of motion 

T% = Q, (6) 

are consequences of the field equations. 



In Brans-Dicke gravity, the modification was introduced due to the lack of compatibility of Einstein's theory 
with Mach's principle. But in Sdez and Ballester [3] a scalar-tensor theory of gravity was introduced in which 
metric is coupled to scalar field. Here the strength of the coupling between gravity and the field was governed 
by a parameter uj. With this modification, they were able to solve a 'missing-mass problem'. Several aspects 
of Sdez & Ballester scalar-tensor theory of gravitation in relation to Bianchi cosmological models have been 
explored in literature. 



The energy-momentum tensor Ty for a perfect fluid has the form 

Ti] = (p + p)u^UJ - pg^j , (7) 

where p is the thermodynamical pressure, p the energy density, Ui the four- velocity of the fluid satisfying 

g.ju'u^ = 1. (8) 

In co-moving system of coordinates, we have Ui = (0, 0, 0, 1). For the energy momentum tensor ([7]) and Bianchi 
type-/ space-time ([1]), Einstein's modifled field equations ([5]) yield the following five independent equations 

A B AB „ 1 ^fe;2 

^ + ^ + ^ = -8.p+-^V, (9) 
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! + § + = 

B C BC „ 1 ,i, :9 , \ 

5 + C + BC=-«'^^+2-^^' (1^) 

is AC BC ^ 1 ^^^-2 

•• ■ ( A B C\ TO(/)2 ^ ^ 

'^^-^(l + ^ + cj+^ = «' (1^) 

where an over dot denotes derivative with respect to cosmic time t. The law of energy-conservation equation 
{t!^ = 0) gives 

P+(p + .)(^ + f + g)=0. (14) 

It is worth noting here that our approach suffers from a lack of Lagrangian approach. There is no known way 
to present a consistent Lagrangian model satisfying the necessary conditions discussed in this paper. 

The spatial volume for the model ([1]) is given by 

y3 ^ ABC. (15) 

We define a = {ABC)^^ as the average scale factor so that the Hubble's parameter is anisotropic and may be 
defined as 



a 






B 


C 


a 








'c 



H-Z = -o\- + ^ + -\- (16) 
We also have 

H=^{H,+Hy + H,), (17) 

where ^ ^, Hy ^ ^ and = %. 

The deceleration parameter q, the scalar expansion ^, shear scalar and the average anisotropy parameter 
Am are defined by 

ad a , ^ 

1=-— = 772' 18 



(20) 



i=l 

where AiJ, = Hi — H(i = x, y, z). 



3 Solution of the Field Equations 

The field equations ([5|)- P^ are five equations necessitating six unknowns A, _B, C, p, p and 0. One additional 
constraint relating these parameters is required to obtain explicit solutions of the system. We consider the 
deceleration parameter as time dependent to get the deterministic solution. 



We have revisited the solutions of Kumar and Singh f25^. Now subtracting ([9]) from (fTO|) . ^ from (jlip . ([T0| 
from (jlip and taking second integral of each expression, we obtain the following three relations, respectively 

^ = di exp (xi j a-^dt\ , (22) 
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^ = rfa cxp (^X2 J a-^d?j , (23) 

^ = exp (^X3 J a-^d?j , (24) 
where di, ^2, d^, xi, X2 and x^ are integrating constants. 

From above Eqs. (|22 p - ((24|) . the metric functions A{t), B{t) and C{t) are expUcitly obtained as 

(25) 



v4(<) = flioexp j a ^dt^, 
B{t) — a2aexp ^62 j cT^dt 



(26) 

C{t) = 030 exp ^63 j a-^d?j, (27) 



where 



oi = (^1^2)3, 02 = (—j , a3 = (d2d3)"3, 

^1 = ^(a;i + a;2), 62 = ■^(a;3 - a;i), 63 = ~^(a;2 + 2:3). (28) 
It is worth mention here that these constants also satisfy the fohowing two conditions 

010203 = 1, 61+62+63 = 0. (29) 

The second integral of (fT3| leads to 

'h{k + 2) 



a^^dt 



(30) 



where h is a constant due to first integral while the constant of second integral is taken as zero for simplicity 
without any loss of generality. 

3.1 Case 1: when a{t) = 

Following Saha et al. [40] and Pradhan and Amirhashchi [41], we take following ansatz for the scale factor, 
where increase in term of time evolution is 

a{t) = VFe^, (31) 

where n is a positive constant. This ansatz generalized the one proposed by Amirhashchi et al. [42]. The 
motivation to choose such scale factor ([31]) yields a time-dependent DP. Recently Yadav [43] has studied such 
type of scale factor in Bianchi-V string cosmological model at late time acceleration. 

Substituting Eq. ^ into Eq. ([H]), we find 

9=r^-l- (32) 

[n + i)^ 

From Eq. we observe that g > for t < \/2n — n and g < for t > \/2n — n. It is observed that for 

< n < 2, our model is evolving from deceleration phase to acceleration phase. Also, recent observations of 
SNe la, expose that the present universe is accelerating and the value of DP lies to some place in the range 
— 1 < g < 0. It follows that in our derived model, one can choose the value of DP consistent with the observation. 
Figure 1 graphs the deceleration parameter [q) versus time which gives the behaviour of q from decelerating to 
accelerating phase for different values of n which is consistent with recent observations of Type la supernovae 
(Riess et al. [26]; Perlmutter et al. [27]; Tonry et al. [31]; Riess et al. [32]; Clocchiatti et al. [33]). 
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> t 



Figure 1: The plot of deceleration parameter q versus t 



Using Eq. (|3T|) in (|25 |) - ([27| . we get the following expressions for metric coefficients 

A = ai{Vv^e*-) exp 



B = 02 (Vt^) exp 

C = 03(^^6*) exp 
Again substituting ([5T|) in ((5n|) . the scalar field is obtained as 

'h{k + 2) 





f {f'eY 


3 , 




f (t"e*)- 


3 , 




f (i"e*)- 


3 , 



Eq. ((55)) implies that 



(33) 

(34) 
(35) 

(36) 
(37) 



Substituting ([55)) - ([57)) in ([TT|) and (|12p . the expressions for thermodynamical pressure (p) and energy density 
(p) for the model arc obtained as 



n 3 /n 
i2 " 4 



uh-" - /33 (re') 



j.n A \ —3 



J(^l)^Q^/^^+/^.)(^v)-3, 



(38) 
(39) 



where 



f3i = bl + bl + bl 

P2 = bib2 + ^2^3 + ^3^1, 
/53 = &2 + &3 + ^2fo3- 



(40) 
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n = 0.5 — - — n = \ n=2 



Figure 2: The plot of energy density p versus t. Here uj — \, h — = 



In view of (|29|) . it is observed that the solutions given by (p3| -(j40 |) satisfy the energy conservation equation 
([T4|) identically and hence represent exact solutions of the Einstein's modified field equations (|9|)- (fT3| . From 
Eq. ([5^ . it is observed that p is a positive decreasing function of time and it approaches to zero as f — > oo. 
This behaviour is clearly depicted in Figure 2 as a representative case with appropriate choice of constants of 
integrations and other physical parameters using reasonably well known situations. 

The rate of expansion Hi in the direction of x, y, and z read as 

H, = \{-^+l)+b,{eeYK (41) 

ffy = ^(^ + i)+fe2(tV)-i, (42) 



i?.= 2(7 + 1; +^3(i'V)-5. (43) 
The Hubble parameter, expansion scalar and shear of the model are, respectively given by 



(45) 

^2 = i/3i(i"e*)-3, (46) 

where /3i is given in Eq. ([Ml. 

The spatial volume iV) and anisotropy parameter {A„i) are found to be 

V^ = (t"e*)i, (47) 
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Figure 3: The plot of anisotropic parameter Am versus t. Here /?i = 1 



A™ = -/3i(tV)-3(- + l) . (48) 

From Eqs. (|47)) and (|45|) we observe that the spatial volume is zero aX t — Q and the expansion scalar is infinite, 
which show that the universe starts evolving with zero volume at t = which is big bang scenario. From Eqs. 
([55)1 - ([55)1 ■ we observe that the spatial scale factors are zero at the initial epoch t = and hence the model has 
a point type singularity |46) ) . We observe that proper volume increases with time. 

The dynamics of the mean anisotropic parameter depends on the constant Pi — b\ + h\-\'h\. From Eq. p8|) . 
we observe that at late time when t — cx), Am 0. Thus, our model has transition from initial anisotropy to 
isotropy at present epoch which is in good harmony with current observations. Figure 3 depicts the variation of 
anisotropic parameter {Am) versus cosmic time t. From the figure, we observe that Am decreases with time and 
tends to zero as i ~> cx). Thus, the observed isotropy of the universe can be achieved in our model at present 
epoch. 

It is important to note here that limj^o (^) spread out to be constant. Therefore the model of the universe 
goes up homogeneity and matter is dynamically negligible near the origin. This is in good agreement with the 
result already given by Collins [44) . 




3.2 Case 2: when a{t) = Vte* 

Following Amirhashchi et al. [42], we consider the following ansatz for the scale factor, where the increase in 
terms of time evolution is 

a{t) = \/te*. (49) 
By the above choice of scale factor yields a time dependent deceleration parameter. Using ([49]) into ([TSl) . we find 

'^-(TTIF^'- ^''^ 
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Using Eq. (HHl) in we get the following expressions for metric coefficients 



A = ai{Vte^) exp 

B ^ a2 (Vte* ) exp 

C = a3{Vte^) exp 
where the constants are already defined in previous section. 





f iter 


^dt 




f iter 






f iter 


2dt 



(51) 

(52) 
(53) 




(0 = CO = 0.5 



(0=1 



Figure 4: The plot of energy density p versus t. Here h = P2 — ^ 



Again substituting (|^ in ((5D|) . the scalar field is obtained as 

'hik + 2) 



Eq. ([54]) implies that 



(te*)~5dt 



(54) 
(55) 



Substituting ([ST])- ([551) in (HB) and ((T2|), the expressions for thermodynamical pressure (p) and energy density 
ip) for the model are obtained as 



P= jQ + l) +(lu;h' + P2]iter' 



(56) 
(57) 



where ^2 and P3 are already given in previous section. 
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Figure 5: The plot of anisotropic parameter Am versus t 



2 

1 ^l 



In view of ((29)) . it is observed that the solutions given by ((5T|) - (|57p satisfy the energy conservation equation 
identically and hence represent exact solutions of the Einstein's modified field equations (P|)- p^ . From Eq. 
(|57p. it is observed that p is a positive decreasing function of time and it approaches to zero as t — > oo. Figure 
4 shows this behaviour of p. 

The rate of expansion Hi in the direction of x, y, and z read as 

iJ, - i fi + l^ +6i(te*)-i, (58) 

i?, = - + +fe2(te*)-i, (59) 

ff.-^^ + l)+&3(te*)-i. (60) 
The Hubble parameter, expansion scalar and shear of the model are, respectively given by 

H^Q + l). (61) 

. 3 /I 



^\p^{te')-\ (63) 
The spatial volume {V) and anisotropy parameter (A^) are found to be 

V^{te')i, (64) 

An = ^/3i(ie*)-3Q + l') . (65) 
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where /3i is already defined in previous section. 

This model has the same properties as described in Casel. This is a particular case of previous one. 
3.3 Case 3: when a{t) = -1 + 




Figure 6: The plot of deceleration parameter q versus t 



In this case, we consider the following ansatz for the scale factor, where the increase in terms of time evolution 



a{t) = -]+t\ 



which yields a time dependent deceleration parameter. Using (j66p into (jlSp . we find 



V 2t3 + 1 



(66) 



(67) 



From Eq. (|S7|. we observe that q = ~2 a.t t = and g = — O.Sasi— ?>oo. Figure 6 depicts the variation of 
deceleration parameter versus cosmic time t. From figure we observe that q = —2 at t = and increases very 
rapidly to reach its maximum near zero and then decreases with time and ultimately it becomes constant (—0.5) 
as t — )■ oo. The values of the deceleration parameter separate decelerating {q > 0) from accelerating {q < 0) 
periods in the evolution of the universe. Determination of the deceleration parameter from the count magnitude 
relation for galaxies is a difficult task due to the evolutionary effects. The present value qo of the deceleration 
parameter obtained from observations are —1.27 < qo <2 (Schuecker et al. [15]). Studies of galaxy counts from 
redshift surveys provide a value of go = O-li with an upper limit of go < 0-75 [45]. Recent observations show 
that the deceleration parameter of the universe is in the range —1 < g < i.e go ~ —0.77 [46]. Thus we see 
that the value of q at present epoch in our derived model is very near to the value obtained by recent observation. 



Using Eq. 



in (|55)) - ([55]) . we get the following expressions for metric coefficients 



A = ai ( -- + <^ ) exp 



6i ( ) dt 
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Figure 7: The plot of energy density p versus t. Here h = P2 — ^ 



B = 02 ( -- + r ) exp 



C = a3 ( -i+t" )exp 



I dt 



t'\dt 



where the constants are already defined in previous section 3.1. 
Again substituting in ([50]) . the scalar field is obtained as 



h{k + 2) 



+ r 



dt 



Eq. dni) imphes that 



t 



+ r 



(69) 
(70) 



(71) 



(72) 



Substituting (|68l) - ([72|) in ([TT]) and (fT2)) . the expressions for isotropic pressure (p) and energy density (p) for the 
model are obtained as 



p = [ ^ivh^ - 13, 



1 



where /32 and ^3 are already given in section 3.1. 



1 



+ r 



3(2^ + 1)^ 

i2(t3-l)2' 

3(2^2 + 1)2 

t2(i3 - 1)2 ' 



(73) 
(74) 



In view of (P^)) . it is observed that the solutions given by (|55)) - (|74p satisfy the energy conservation equation 
([H]) identically and hence represent exact solutions of the Einstein's modified field equations ©-([ISl)- From Eq. 
(j74p . it is observed that p is a positive decreasing function of time and it approaches to zero as t — s- cx). Figure 
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Figure 8: The plot of anisotropic parameter Am versus t 
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7 depicts the variation of energy density versus cosmic time t. We observe from the figure that p is a positive 
decreasing function of time and it approaches to zero at late time. 

The rate of expansion Hi in the direction of x, y, and z read as 

Hy= \Z +bA-\ + t') ' (76) 

The Hubble parameter, expansion scalar and shear of the model are, respectively given by 

_ (2t3 + l) 

_ 3(2^3 + 1) 

.^4/3r(-i + .^)" (80) 
The spatial volume (V) and anisotropy parameter (A™) are found to be 

3 



V=[-l+t'\ . (81) 
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where /3i is already defined in previous section. 

It is observed that the spatial volume is zero at t — 1 and the expansion scalar is infinite, which shows that 
the universe starts evolving with zero volume at initial epoch t — 1 with initial rate of expansion. The scale 
factors also vanish at t = 1 and hence the model has a Point type singularity [17] at the initial epoch. The 
pressure, energy density, Hubble's factors and shear scalar diverge at the initial singularity. Figure 8 is the 
plot of the anisotropic parameter (Am) versus cosmic time t. From the figure, we observe that Am decreases 
with time and tends to zero as i — J' oo for all values of /3i. Thus, the observed isotropy of the universe can be 
achieved in our derived model at present epoch. The shear tensor also tends to zero in this model. Thus the 
model represents shearing, non-rotating and expanding model of the universe with a big bang start approaching 
to isotropy at late times. 



4 Concluding Remarks 

In this paper we have studied a spatially homogeneous and anisotropic Bianchi type-I space-time within the 
framework of the scalar-tensor theory of gravitation proposed by Sdez and Ballester [3]. The field equations 
have been solved exactly with suitable physical assumption. The solutions in all three cases satisfy the energy 
conservation Eq. p4)) identically. Therefore, exact and physically viable Bianchi type-I models have been ob- 
tained. To find the deterministic solution, we have considered three types of scale factors which yield three 
different time dependent deceleration parameters. 

In Case 1, we choosed the scale factor a{t) = -\/t"e* which yields a time dependent deceleration parameter 
9 = {n^t)'^ ~ Since all the scale factors vanish at initial moment, hence the model has a big bang singularity 
at t — 0. This is a Point type singularity. This case represents a model which generates a transition of the 
universe from early decelerated phase to the recent accelerating phase which is in good agreement with the 
recent observations (Caldwell et al. [15]) (see. Fig. 1) whereas in Pradhan et al. [33] only the evolution takes 
place in an accelerating phase. The model has also transition from initial anisotropy to isotropy at late time 
(i.e. present epoch) which is in good harmony with current observations (see. Fig. 3). For different choice of n, 
we can generate a class of cosmological models in Bianchi type-I space-time. It is observed that such models are 
also in good harmony with current observations [26, 27, 31-33]. 

In Case 2, we considered the scale factor a(t) = Vte^ which yields a time dependent deceleration parameter 
Q — (i+ty^ ~ 1- This is a particular case of previous one. This model has the same properties as obtained in 
Case 1. It is worth mention here that the anisotropic parameter becomes negligible at i — ^ cx) which, in turn, 
implies that the derived model isotropizes at present epoch. 

In Case 3, we considered the scale factor a{t) = — \+t^ which yields a time dependent deceleration parameter 

q = —2 ^ 2t^+i ) • ^^^^ '"^^^ accelerating universe (see. Fig. 6) at present epoch. In this case, the 

anisotropic behaviour of the universe dies out on later times and the observed isotropy of the universe can be 
achieved in derived model at present epoch (see. Fig. 8). The model has a point type singularity at initial epoch 
t = 1. 

It is worth mentioned here that in all above three cases, their anisotropic solutions may enter to some 
isotropic inflationary era. All the above three models represent expanding, shearing and non-rotating universe, 
which approach to isotropy for large value of t. This is consistent with the behaviour of the present universe as 
already discussed in introduction. If we set fc = 0, then the solutions reduce to the solution in general relativity 
and the Sdez-Ballester scalar-tensor theory of gravitation tends to standard general theory of relativity in every 
respect. In literature we can get the solutions of the field equations of Sdez-Ballester scalar-tensor theory of 
gravitation by using a constant deceleration parameter. So the solutions presented in this paper are new and 
different from other author's solutions. 

Finally, our solutions may be useful for better understanding of the evolution of the universe in Bianchi-I 
space-time within the framework of Sdez-Ballester scalar-tensor theory of gravitation. The solutions presented 
here can be one of the potential candidates to describe the observed universe. 
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